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^ Abstract: Born-lnfeld nonlinear electrodynamics with point singularities having 

■ both electric and magnetic charges are considered. Problem of interaction between 

. the associated soliton dyon solutions is investigated. For the case of long-range inter- 

O ■ action at first order by a small field of distant solitons we obtain that the generalized 

Lorentz force is acted on a dyon under consideration. Short-range interaction between 
0^ , two dyons having identical electric and opposite magnetic charges is investigated for 

. an initial approximation. We consider the case when the velocities of the dyons have 

equal modules and opposite directions on a common line. It is shown that the associ- 
Oh. ated field configuration has a constant full angular momentum which is independent 

^ . of the interdyonic distance and their speed. This property permits a consideration 

of this bidyon configuration as an electromagnetic model of charged particle with 
^ . spin. We numerically investigate movement of the dyons in this configuration for 

the case when the full electric charge equals the electron charge and the full angu- 
lar momentum equals the electron spin. It is shown that for this case the absolute 
value of relation between electric and magnetic charges of the dyons equals the fine 
structure constant. The calculation gives that the bidyon may behave as nonlinear 
oscillator. Associated dependence of frequency on the full energy is obtained for the 
initial approximation. In the limits of the electrodynamic model we obtain that the 
quick-oscillating wave packet may behave like massive gravitating particle when it 
move in high background field. We discuss the possible electrodynamic world with 
the oscillating bidyons as particles. 
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1. Introduction 

Solitons in nonlinear electrodynamics may behave like real particles. This analogy 
with classical particles is manifested when we consider a long-range interaction of 
solitons by a perturbation method 0. In this case a soliton is subjected to the 
Lorentz force in the first order by a field of distant solitons. The second order may 
give the soliton's trajectory in the form of geodesic line for some effective Riemann 
space with metric depending on the field of distant solitons. Thus we have also 
an analogy with gravitating particles. Moreover, light beams distortion under the 
action of some given field may appear as gravitational distortion [Q. These properties 
(in particular) provoke interest to nonlinear electrodynamic models in the context 
of possible unifying description for the real material world. In the framework of 
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this approach it may hope to solve the problem of unification electromagnetism and 
gravitation. 

Here we shall consider an interaction of singular solitons in Born-Infeld nonlin- 
ear electrodynamics^ . The investigation of this nonlinear electrodynamic model with 
point singularities of field is presented in the article , where we have considered the 
singularities with electric charge. But the model's system of equations for electro- 
magnetic field has an exact solution which has non-zero divergence for both fields of 
electric D and magnetic B inductions at the singular point. That is, this point for 
the solution has both electric and magnetic charges. The particle with both electric 
and magnetic charges had been named by J. Schwinger as dyon (see [0). Because the 
singular solution looks like point charged particle we name this field configuration 
as dyon. In the present article we consider an interaction between these dyons and 
we consider some quick-oscillating two-dyons field configuration. We also consider 
like-gravitational interaction of a quick-oscillating wave packet with given field. 



2. Basic relation for field 

Let us state the basic relation for the Born-Infeld nonlinear electrodynamics with 
singularities [^] generalized to the case with both electric and magnetic charges. We 
shall use an inertial coordinate system in Minkowskian space: {x^} (the Greek 
indexes take values 0, 1, 2, 3 ). That is, components of metric is independent of time 
and \goo\ = 1, goi = (the Latin indexes take values 1,2,3). In this case (see 
also 0) we can write the following nonlinear Maxwell system of equations: 



divB 
divD 
(9o B + rotE = 
^ doB- rotH = 



= Anf 
-- Anf 
-An J 
-An J 
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^It will be recalled that M. Bom and L. Infold had been considered in their article the 
electrodynamic model which follows from an action that had been proposed by A.S. Eddington in 
his book g. 
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and are components of electric and magnetic singular currents such that 

N N 

V n=l V l^'l „=i 

iV iV ^^-^^ 

V n=l V 1^1 „=i 

n n 

d and 6 are electric and magnetic charges of ra-th singularity, 
a = a(x ) IS a trajectory of it and V = -r—^ ■ 

yXJb 

Here we use the definition for the three-dimensional 5 -function which is suitable for 
discontinuous functions /(x) : 

j /(x) 5(x - a) {dxf ^ hm ^ I /(x) da ^ ^/(x))^ , |a| = | da . (2.4) 

n I I I n I ^ 71 

o \ a J a 



where f2 is a region of three-dimensional space including the point x = a , 
a is a closed surface enclosing this point, da is an area element of the surface a , 
I a I is an area of the whole surface a . 
Let us define the following functions: 

'H = C + a^Ti-^ = Z + a^B-ll = a'^T^^ + l , 

P = (D X B) , -pt ^rj^Qi ^ { ■ ) 

where T^^ are components of symmetrical energy-momentum tensor 
(see where we designated it as T^" ). 
Then, from relations (|2.2|) we can obtain the following relation (see also |^): 

;J-(D-a2pxB) = [(1 + «2 J) _ ^2 (D.B)B] 

, (2.6) 

l-{B + a^Vx-D) = l[(l + a2D2)B_^2 (d . B) D] 
rt ri 
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OH 
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dB 




n 



where H = Jl + a^B^ + B'^) + V 



;i + a2D2)(l + a2B2) -a4 (D.B)2 . (2.7) 



Using relations ( p. 61) we can consider system (|2.1|) as the system of equations 
for fields D, B . This representation is best suitable for investigation of the singular 
dyon solutions. 

Let us introduce the two electromagnetic potentials , . In our case the 
potentials have singular line for each singular charged point (in three-dimensional 
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space). We must exclude such line from space when consider any ordinary differential 
field model. The alternative way is connected with using distributions or generalized 
functions. As described by P. A.M. Dirac for monopoles |]10[, in this case we must 
include some distributions into definitions of the potentials through derivatives. But 
here we shall adhere to the first way associated with exclusion of the singular lines. 
Note that the singular currents into equations ( p.lD set boundary conditions at the 
singular points (see 0). In contrast, we shall take natural boundary conditions at 
the singular lines outside of the singular points. 

Thus we define the potentials outside of any singular set with help the following 
formulas: 



where £0123 = \9 
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We have the following definition in three-dimensional designations: 



[2.9) 
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-VAo + do A 
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(2.10) 



(b) 



From the basic equation of the model outside of the singularities 
1 d 



up 



f2.111 



we can easy obtain also the following equation for the potential A^ outside of the 
singular set: 



where 



dxf^ 



, 



(2.12) 



(Jf^uap ^ gpa gUp _ gpp ^u^ _ ^2 ^jrf^u jr-p ^ jpu jap^ ^ ^2.13) 
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3. Boundary conditions at the singular points 



Let us integrate system of equations over a small four-dimensional space region 
including only n-th singular point. Using the partial integration in four- dimensional 
space we obtain the following conditions: 



J J ff"" da^ -ATxjfdndx^ = Q, j j T^"" da^ -AttJj'' dn 



dx^ 







(3.1) 



where is a three-dimensional region including only n-th singular point, 

da is external directed element of the closed (two-dimensional) surface 
a enclosing the singular point (see also [§]), 
dQ = J\g\ (dxY , dao = - V ■ do- , 



n 1 

do"^ dx° — ~ Q ^fiiyap dx" A dx'^ A dx^ are components of four- vector. 

The surface a is rigidly coupled with n -th singular point and move with it together. 
Relations ( p.l| ) will be satisfied if we have in three-dimensional designations 
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(a) 



(3.2) 



(b) 



Because the surface a may be arbitrarily small, relations (|3.2| ) are boundary condi- 
tions at the n-th singular point. 

There is the attractive idea, proposed by A.S. Eddington 0, about an invariance 
of theory under permutation of points of space. We can apply this idea to the singular 
points. Let the model be invariant under permutation of any two singular points in 

ni 712 ni n2 

the sense of change their charges d^ d , b^ b- Because the theory is invariant 
under change of charge's sign, in this case we have for any singular point 

71 — 71 — 

d= ±d , b= ±b , (3.3) 

where d and b are some positive constants of space. 

Thus we have the theory with three dimensional constants a , d and b . For 
the suitable dimensional system we can take a = 1 and d = 1 or 6 = 1. Thus we 
have the relation d/b as the single dimensionless constant of the theory. Below (in 
section |9|) we shall connect this relation with the fine structure constant. 
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4. Dyon solutions and singularities 



System of equations with relations (|2.6| ) has the following exact dyon solution 
in a cartesian coordinate system {y^} : 



where 



±6 



.0 - dMr) 
io = -bMr) 
d = ± d 



y^ b 

A = 
A = 



b 



(b) 



(4.1) 



&0(y) 



yvi 



dr' 



,/4 



a 



{d' + F 



(4.2) 



The vector function (y) may be of two types: with infinite singular line and with 
semi-infinite one. Its components have the following forms: 

1 y'^y^ I y^y^ 



or 



r 

1 (r + y^) y"^ 
r 



r 



,03 = 

1 (r + y^) y^ 
r p^ 



(4.3) 







(b) 



where p = \l {y^Y + {y'^Y 

Function ( ^4.3| a) has the infinite singular line coinciding with the axis y^ and 
functions (Ob) has the semi- infinite singular line in positive direction of y^ . In 
a spherical coordinate system {r, if} the both functions </> (y) have non-zero 
(^-components only that can be written in the forms 

(a) or (6) 

cot 7? ^ I d (4-4) 

005 = 0(p = cot - 

r r 2 

accordingly. 

At the singular line the potentials A, A are devoid of defined direction and their 
absolute values are infinity. At the singular point r = the vectors E, H, D, B are 
devoid of defined direction and absolute values of D, B are infinity. 



The full electromagnetic energy (see (|6.3| ) below) of dyon solution (|4.1|) is 



(4.5) 



The dyon's energy has a space locahzation region. Let us denote the sphere- 
enclosed dyon's energy as £' for radius of the sphere r' . Center of the sphere is at 
the origin of the coordinates {?/*} . The numerical calculation gives that £' = 0.5 £ 
for r' = f and £' = 0.95 £ for r' = 10 f. 

With help of Lorentz transformation, shift, and rotation of the coordinate system 
{?/*} we can obtain the following moving dyon solution in the coordinates {x^} : 

= (L] - e'i^ Vj Hi) j Vl - V2 

f-. . .., X / , , (a) 

3' = [L] & + e*^' Vj Rij I Vl - V2 

5^ = (p. - e'^^ Vj Bj) I VI -V2 

V = (Z} W + e'^^ Vj Di) I VI- V2 

where D! , R\ R\ K\ Au, Au are defined by formulas 
with = Dj {x^ — V^ — al) and = TV-j , 
L^^ is the Lorentz transformation matrix, U- Lj = 6] , 

^.0 - 



(b) (4.6) 




(4.7) 



'1 - V2 

Oq are components of an initial position, TZ^j is a rotation matrix. 
5. Variational principle for two potentials 

Here we propose some action such that the associated variational principle gives 
system of equations ( |2.1| ) with relations ( pT6| ) and boundary conditions (|3.2|). This 
action has the following form: 



S 



J[{C + 1) +27ra'{fA,-rA^)]^{dx)^ , (5.1) 



«2 



where £=-7^+y(E-D + H-B) , (5.2) 



H = H(D, B) according to 
E, B are represented by A^ and D, H are represented by A^ with ( |2.10| ). 



Using definition (|2.5|) we can easy write also the following expression for C : 

1 
2 



£ = -i(£ + £) . (5.3) 
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Thus the action S_ is invariant under general transformation of coordinates. 

If we substitute expression ( ^.31) for £ into the action S_ and take into account 
that £ = C{d^Ai,), C = C^d^A,^) , then independent formal variation of eight com- 
ponents of the potentials A^, for the principle 6S_ = gives the following eight 
equations: 



d 



\g\ ^x^' 



dC 



d{d,A^ 



47r/ 



d 



\g\ dx^" 



dC 



d{d,A_ 



(5.4) 



For simplicity we didn't consider here the singular lines of the potentials. 

Formally we can consider system of equations (|2.1|) with relations (|2.2| ) as the 
system of eight equations for the eight unknown functions A^, A^, ( (|2.1| aO, (|2.1| y) for 
Afj^ and (|2.1| a^0,(prT^") for A^). Such system of equations agree with system ( ^.4|) . 

Now let us consider the action S_ with the function C in form (|5.2| ). In 

this case the dependence of Lagrangian on the derivatives of potentials differs from 
one that we have for the function C in form ( ^.3|) . Thus the appropriate variational 
principles with the function C in forms ( |5.2| ) and ( |5.3| ) are different. 

Substituting definitions ( |2.10[ ), (|2.3|), ( |2.4| ) into ( |5.1| ) we obtain the following 
expression for the action that doesn't contain 5 -functions: 



(5.5) 



blAo + V-A])dx" 



S = J (£ + l)v^(dx)^ + 



a 



VAo - doA) • (V X a) + ((9oA - VAq) • (v x A 



(5.6) 



n 

a 




Note, though potentials A^ , A^ (|4.6| c) 
for the dyon solutions are infinite at the sin- 
gular lines, its averaging at the singular points 
{A^)n , {A^)n are finite. Thus the action S_ is 
finite for the dyon solution. 

Because definition of the potentials ( p.lO| ) 
may be used only outside of the singular set, 
expression for C ( p.6[ ) is true only outside of 
the singular lines. Let us enclose all singular 
set in a multi-tuple connected surface a with 
external (relative to the singularity) surface 
element da- . Let this surface be composed 
of parts of surfaces a for the segregated 

singular points and multi-tuply connected tubular surface a (in three-dimensional 
space) for parts of the singular lines outside neighbourhoods of the points. For 



dcr = dcr 

Figure 1: Part of the surface a for 
the dyon field configuration with semi- 
infinite singular line. 
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example, the surface a for the semi-infinite singular line can have image as it is 
shown in figure |l|. 

Let us designate the action S_ in space outside of the singular set as S_ . We can 
obtain variation of this action with partial integration. Thus we have 



ss 




do'D - V X 



dB 



■6A 



doB + V X 



m-c 



6A 



a 
~2 



+ 



+ 



D ■ d&) SAo - (B • do-) SAo\ + 271 a^Y. [d {^^o}^ - b {SAo) ^ 
— X dCT-y (Ex dCT + Bdao) ■ 6A-2na^Y.b^ ■ {^^ 



X dcr — — (H X dcr — D dcrn) 
OB 2 ^ ^' 



■5A + 2Txa^YA^ ■ {^^ 



(5.7) 



+ 2-na^Y. 



d E + V X B +6 H- V X D 



n 

■ oa 



£ dcr ■ (5a , 



where Vl is three-dimensional space without the region bounded by surface a, 

da 



dcTo = — V • dcr 



V 



dxO 



a is a points of the surface a 



Now we contract the surface a to the singular lines outside neighbourhoods of 
the N singular points. Next we contract the surfaces a to the N singular points. 
As result we have 5S_^ 5S_. 

Let us find stationary conditions for the action S_. By the general principle of 
calculus of variations [|ll], we can take (5a = (5a = at first. We can also take 
continuous variations (5A^, (5^1^ and make, firstly, 5Afj_ = 0, (5^4^ = at the singular 
lines. Then, according to the first line in expression (|5.7| ), we have 



(9oB + V X ^ = 



dD 



a'do^-V X 1^ = 



(5. 



in the space outside of the singular lines. Using definition of the fields B, D, E, H 
through potentials (|2.1CI|) , we obtain from equations (|5.8|) that 



E=l^ 
a2 dB ' 



a2 OB 



(5.9) 



Now we take 6Afj_ ^ and (5^4^ 7^ at the singular lines for continuous variations 
(5^4^, (5y4^ . We assume that the fields D, B, E, H are continuous at the singular lines 
outside the segregated singular points. Then we obtain boundary conditions ([37 
from expression ( pT7| ) (line 2-4), using formulas (|5.9|). 
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Finally, we take 5a 7^ , 5a^ . It would appear reasonable that at the singular 
lines 

\im fc do- =0 . (5.10) 



Then, from the last line in (|5.7|) we have the following condition at the n-th singular 
point: 



(i(^E + V X Bj +6(^H- V X Dj^ =0 . (5.11) 

n 

This condition generalize one for purely electrical singularity {b = ) , which had 
obtained in the article 0, to the case of dyon singular point. Using definition of 
singular currents ( |2.3| ) we obtain also the following form of the condition: 

F,.f-^e,,^,rr = . (5.12) 



Thus the variational principle for action S_ (|5.1| ), (|5.5|) gives system of equations 
[1]) and point boundary conditions (|3.2|) , ( |5.11|) . 



Notice that here we don't use a concrete form for the function 7i(D, B) . Thus 
this derivation is suitable for any electrodynamic models. 

6. Conservation laws and comments about dimensions 



Using equations ( p.l|) with formulas ( |2.6|) and condition (|5.12|) we can check directly 



the following conservation law for symmetrical energy-momentum tensor (in cartesian 
coordinates) : 

where 

TOO = «-2 {n-\) , 

yoi ^ e'jiDjBi = e'^^EjHi , (6.2) 
= [D ■ E + B ■ H - T°°] - {D' + B' W) . 

Let us introduce the following notations for full energy and momentum in the 
region VL of three-dimensional space VL ■ 



iy^^^^ ' (6.3) 



where = or P = D x B . 
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Let the region Q may be moving, that is we have Q = Q{x^) . Then, by 
integration of conservation law (|6.1|) over the region Q we obtain the following 
integral conservation law for energy-momentum: 

^ lV-d&- I in -I) dan , (a) 



dxO 



(6.4) 



where a is the surface that bounds the region Q , 
ddj is the external element of this surface, 
ddo = — V • dcr , V is the velocity of surface's point. 

It is evident, from ( |6.1| ) we have also the following conservation law for angular 
momentum tensor: 



, (6.5) 



where 



M^^^ = T'^^ xf - Tf'P x" . (6.6) 



We introduce the following notation for vector of full angular momentum in the 
region Q : 

M = {V X x) df] . (6.7) 

Let us take up somewhat the dimensions of the quantities. Because we don't 
explicitly use the speed of light c , the temporary coordinate x*^ has a dimension 
identical to that for the space coordinates , that is dimension of length. (But we 
can make substitution x^ ^ ct for all time.) For the same reason the momentum 
has a dimension identical to that for the energy. And dimension of the full angular 
momentum is [M] = [P] [x] = [S] [x] . In this case the dimension of mass is identical 
to that for the energy. 

7. Method for investigation of interaction between dyons 

We shall consider system ( |2.1| ) with relations ( |2.6| ) as system of equations for the 
fields D, B . Let us write this system in the following formal operational form: 

AfV = , where ^= (g) • (^-l) 
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We can search a solution of equation ( |7. 1|) in the form of the following formal 
series: 

k 



(7.2) 



1=1 



where P*-"-* is an initial approximation, 
V'^^ is fc-th correction, 
P*^'^) is fc-th approximation. 
Let us expand the operator into equation ( [7.1|) near to the k -th approxima- 
tion P'^'^) and keep only linear terms by the correction T>^^^^ . Thus we obtain the 
following equation for the correction: 



(7.3) 



where M' is Frechet's derivative for the operator M . 

Equation (|7.3|) gives the following iterative process, which is called Newton's 
method for operators [p!2| : 

= _ jTv-' [r'W]}"^ [UV^^'^] . (7.4) 

An alternative possible way is based on using the principle of contractive map- 
pings. Let us rewrite system ( [7. 1| ) in the following form: 



MV = MV 



where MV = MV - AfV 



(7.5) 



and TVI is the operator represented by left parts of linear Maxwell equations for 
D = E, B = H. 

In this case we can easy obtain the appropriate inverse operator Ai^^ . As result 
we have the following representation for equation (|7. 1|) : 



(7.6) 



If the operator Ai ^ JV is contractive (for some given class of functions including 
the initial approximation), then the iterative process 



(7.7) 



converges to the solution. 



The convergence of iterative processes ( [7.4|) , (|7.7] ) to some solution is defined 



by conditions that depend on initial approximation [|T^. But we can use these 
procedures, even if it doesn't converge. In this case we may have some number of 
iteration for best approximation to the solution. In any case the initial approximation 
D*^") must be sufficiently close to some exact solution. 
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We can take the sum of some solutions as initial approximation. Let us take 

n 

moving dyons V ( |4.6D as these solutions. Thus we have 

N 

p(o)=^:p . (7.8) 

n 

The field configuration D, B for the exact moving dyon solution has six free 
parameters: three components of the velocity and three components of the initial 
position. We can assume that these parameters for the dyons into initial approxima- 
tion ( [7.8| ) are functions of time. This is standard practice for similar problems. But, 
as it can show, if we take that the initial positions depend on time, then boundary 
conditions (|3.2| b) are broken for the dyon field configuration. Thus we may have only 
three components of velocity as dyon's parameters depending on time. 

We believe that each dyon has an optional trajectory x = a(x°) such that 

n 

n- n. / ■ ri ■\ del " • n ■ n — . — . n 

y' = L]{x^-a^) , ^-v , l} = l;.(v) , l} = l;.(v) , 

(7.9) 

where i/* are included into expressions (|4.1| a), (|4.1| b), 

n 

Dj are included into expressions (^4.6| a). 
The trajectories a(a;°) must be obtained. 

In the same manner as for field we can represent the functions a(a;°) by the 
following series: 

^ = SW + a^) + a2) + ... , ^(^)^S(o)+^SO , (7.10) 

1=1 

where a^^^ is the trajectory of n-th dyon for the initial approximation, 
a'^^ is k-th correction to the trajectory, 
a^'^) is k-th approximation for the trajectory. 

If we have obtained the trajectories in initial approximation a^'^\x^) with help 
some method, then we can obtain the first correction T)^^ as solution to equation 
( [7.3| ) or from formula ( |7.7| ) for A; = . This correction J)^^ is a radiation of the 
dyons accelerated according to their trajectories a(°)(a;°) . But the radiation of the 
dyons modifies their trajectories. We describe these modifications by functions a^^ 
that correct the trajectories a'^'') . Thus we must obtain the corrections a*^^ to the 
trajectories a(*^~^) in any step of iteration. In effect, this is a correction of the initial 
approximation (through the correction of the trajectories) for every iteration step. 

A direct calculation gives that all boundary conditions ( |3.2| ), ( p.ll| ) are satisfied 
for initial field configuration ( [7.8| ) with ( [7. 91 ) and optional trajectories of the dyons. 
This is distinction of this method from one that had used in the article 0, where 
we take the four components of the potential as unknown functions. But using 
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the fields D, B as unknown functions, tlie present metliod is more suitable for 
satisfaction to all point boundary conditions in any iteration step. 

In the present article we shall use integral energy-momentum conservation law 
( |6.4| ) for obtaining trajectories of the interacting dyons. 

Notice also the following. Because near n-th singular point absolute values of 

n n 

the fields D , B tends to infinity, for any finite given external field (of other dyons) 
we can define some neighbourhood of the singular point, where this external field 
may be believed to be small. Thus the consideration of any (small) distances between 
the dyons in initial approximation may be worthwhile. 

8. Interaction of dyon with a small given field 

n 

Let us consider the interaction of moving dyon with a small given field T)^ 

n n n n 

(i = l,...,6), (EfiiD, H;^B)in initial approximation. In particular, this small 
field may represent the field of other dyons that are at sufficiently long distance from 
the dyon under consideration. 

71 

In this case we consider the space region VL that include the dyon under con- 

71 

sideration only. The region VL move with ra-th singularity together. We can take 

n n 

a maximal absolute value of the functions in the region VL as the small order 

n 

parameter e for the problem of finding the first approximation in the region Vt . Here 

n 

we consider the case when the region Vt has more size than one for dyon localization 
region (see section (|)). 

n 

Thus for the region Vt we take 

n 

j)i^)=j) + V . (8.1) 
For obtaining dyon's trajectory equation we use formula (|6.4[b). We consider the 

71 n 

case when the field T) is approximately constant into the region VL . Thus we have 

dp{0) dp d 

=^ ' (8-2) 



dx^ dx° dx° w _ y:, 



where £ is defined into ( [4.5|) . Here we take into account that size of the region VL 
is much more than f ([4.2|) . 

Now we expand the components T*-^ as Taylor series near the n-th dyon's field 
configuration. We restrict the consideration to the terms of no more than the first 
order by e . Thus we have 

n . . 

Ti=r' + ^V'+... , (8.3) 
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where Z = 1, 6 , 

According to formulas (|6T 



n . . 



we have 



n . . 



Qrpij 

dDi 


1 




dBi 


1 


u + 


We have also 






d^n _ a' 




dD^dDi H 




d^n _ a" 




dB^dB^ H 








dD^dBi H 



5; 



6i 6^' 



DP 



DP 



dD' dDi 
dD' dBi 



+ BP 



BP 



dB' dD\ 
dB' dBi , 



(8.4) 



- {B, B, + E, E,)] 
[1 + D2) - {Di Dj + Hi Hj)] 



15) 



Di Bj - 6ij (D ■ B) - Bi Dj 
Using these formulas and formulas ( ^.6] ) we obtain from 

n . . 

dT^ 



Ei Hj] 



that 



dDi 

n . . 

dT^ 
'dB^ 

(I 



-SI E' + 5'^ Eg 



SIE' + Oil, 



.6) 



. n . . . n . n . / /n 

-51 W + 5'^ Hq-5lW + (l /r' 



We substitute ( |8.2| ), ( |8.3| ), ( ^.6| ) into relation (|6.4|b) and make the integration for 
its right part. Because the symmetry properties of element of integration, the last 
integral is zero. As result we obtain the following generalized Lorentz equation for 
the trajectory of n -th dyon: 



V 



dxO 




n n \ 

~ n ~ 

B - V X D 



(8.7) 



1 - V2 

As we see, the energy-momentum method for obtaining the trajectories gives the 
value of mass which is equal to rest energy of the dyon. This result fully conform 
with Einstein's principle for equivalence of mass to energy. 

In the article we considered the potentials as unknown functions and used 
the boundary condition of type ( |5.11| ) for obtaining the trajectory equation. As result 
we had obtained another value of mass in initial approximation {\£) than in the 
present article. Generally speaking, we have certain arbitrariness for taking an initial 
trajectory of the singularity. Because a true trajectory of the singularity is defined 
by an exact solution only, we can obtain corrections to the initial trajectory (in first 
and highest iteration steps) such that it will appear as renormalization of mass. But 
the method of the present article is more suitable, because it gives the satisfaction 
as energy-momentum conservation law ( |6.4| ) as boundary condition ( |5.11| ) for any 
iteration step. 
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9. Bidyon or an electromagnetic model of particle with spin 



Let us consider the interaction of two dyon singularities with identical electric and 

12 1 2 

opposite magnetic charges: d = d, b = —b- At first we use a cylindrical coordinate 
system {z, p, (p} such that dyons are on the axis z . We believe also that the dyons 
velocities are directed with the axis z and the problem is symmetrical about the 
plane z = . This configuration is shown in figure 0, where d = ± d , b = ±6. 



d = d 
b = b 



1 

V 



V 



V = V 



2 

d = d 

2 

b = -( 



Figure 2: Disposition of the two dyons in the cylindrical coordinate system. 



We name this disposition of the singularities with appropriate field configuration 
as bidyon. 

According to our iterative method (see section |^ we take the sum of the fields for 
moving dyons ( [4.6| ) with unknown their trajectory a{x^) as initial approximation. 
That is we have 



V 



(0) 



1 2 

v + v 



(9.1) 



In the cylindrical coordinate system we have the following expressions for the fields 

12 12 

D, D and B, B : 



where 



V 



1 


1 


1 




d 


b 


VI- 




1 

Dp 


1 

Bp 


1 




d 


b 


Vi- 




1 


1 

Bip 


V 




b 




Vi- 





z + a 



da 
d^o 



1 

r 



2 


2 

Bz 


1 




d 


-b 


VI- 




2 

Dp 


2 

Bp 


1 




d 


-b 


VI- 




2 

Dp 


2 

Bp 


V 




b 


~d' ^ 


VI- 


\/2 



z — a 



z' - a' 



(9.2) 



VI -v^ ' VI -v^ 

The electromagnetic potentials don't use in our iterative method. But from the 
initial approximation we can make some suppositions about true field configuration 
of the potentials for the appropriate exact solution. It is reasonable to take that the 
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both dyons into ( |9.1| ) have the singular hnes coinciding with the axis z . That is we 
have 



1 


1 


A, 


A, 


d 


-b 


1 


1 




Aip 


~v 





^ Mr) 



( 2 


2 


Az 


A. 


-d 

< 


-h 


2 


2 


Ai^ 









V 



■Mr) 



0^(2;' - a', p, (p) 



(9.3) 



Here we may take the infinite singular lines with function (p^^ of type (4.4a). 
In this case the singular lines for the potential A particularly cancel each other 
outside of the interval between the dyons. This cancelling is complete as 2; — ± 00 

or a ^ . If we also believe that at the singular points the field V of the bidyon 

1 2 

solution are equal to the field for the dyon solution T> or P , we must believe that 
the exact solution has the potential A with the function (p^ of type ( [4.4| b) at the 



singular points. Thus we can presume that the potential A for the associated exact 
solution has the singular line segment between the dyons only. That is we have a 
singular string segment (for A ) with the dyons at the ends. 

Let us calculate the full angular momentum for field configuration ( |9.1| ). We use 
definition ( |6.7| ) for full three-dimensional space. Because of the symmetry property 
of elements of integration into ( |6.7|) , for this case we have Mp = = , 



M, = -^Jv^p (pdzdpdif) . (9.4) 
We can easy obtain the following expression: 

V, = . (9.5) 

Now we introduce new variables of integration (, ip} that appropriate to the 
bispherical coordinate system: 

, a'sinh^ a' sinC 

P=7i::r7 — z:z7 ■ (9-6) 



cosh ^ — cos C cosh ^ — cos ( 

We have the bispherical element of value 

(pd.'dpd^) = . (9.7) 

(cosh 4 — cosQ 

1 2 

We have also the following expression for r , r : 

1 ^ V2a'exp (e/2) 2 ^ ^2 exp (-^72) 

Vcosh ^ — cos ( ' a/ cosh ^ — cos ( 
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Substituting ( p.6|) , ( p.8| ) into ( |9.5| ) and introducing the variable z = cosC, we 
obtain the following expression: 



-2x1 



= -bd ^^ ' (coshe -zj-^l-z^ 



(9.9) 



And substituting (|^), (|9.7|), ( |9.9| ) into ( |9.4| ) we obtain the following expression for 
z -component of the full angular momentum: 



bd 



OO 1 



— TV — OO 



dz 



(9.10) 



(cosh ^ — ^ 

Making firstly the integration over z_ and ^ in the finite limits [— ^, we obtain 



OO 1 



d^ , 

\ (cosh^-z) 

Thus we have 



dz_ = Jim 

^— >oo 



4 I In - — ) sinh,^ 
\ cosh ^ — 1 / 



(9.11) 



= 2bd 



(9.12) 



As we see, the full angular momentum is independent of the interdyonic distance and 
their speed! It is evident that this result is typical for any electrodynamics model, 
i.e. the form of the function 7i(D, B) is immaterial. 

The property, that the full angular momentum is independent of the internal 
movement parameters, permits a consideration of the bidyon as an electromagnetic 
model for charged particle with spin. The charge of this particle is 2 c? and its spin is 
equal to \Mz\ . We set 2d = e, were e is the absolute value of the electron charge, 

h 

2 



and |M, 



Thus we have 



- eh el 
~ 2 ? ~ 2 ^ 



^ = a 



(9.13) 



where a = e'^/h ~ 0.00729735 is the fine structure constant. 

Now we shall find the trajectory a(x°) in initial approximation. For this purpose 
we use the integral energy conservation law for full three-dimensional space. The full 
energy of initial field configuration ( p.l| ) is depend on the distance between the dyons 
a and on the speed V : £ = S{a, V) . On the trajectory a(x°) the full energy is 
constant. Thus we can obtain the trajectory with help the following formula: 



X 



da 



a(0) 



V{d, S) 



(9.14) 



where V{a, S) is the inverse function to the function S{a, V) . 
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It is evident that we have a family of the trajectories for all possible values of 
the energy. 

But in fact, the energy depend on square of the speed. Thus we can obtain the 
acceleration from the function £^(a, V"^) by the following way: 

dxo ~ da ^ dV^ dxo ~ ^ dxo ~ 2 da/ dV^ ' ^ ' ' 

Let us obtain the function S{a, V^) . We shall use a numerical calculation for 
the energy integral. But at first we obtain some interesting approximate formula for 
the energy in the case of very small distance between the dyons and their speed is 
not very close to the speed of light. We substitute the field T)^^^ into expression 
for H (2.7) and use formulas for bispherical coordinates (9.6) and (9.8). We obtain 
that as a' ^ the main term of the radicand is {'Pipf , which is proportional to 
(a')~^ . Thus for very small values of a' we have the following formula: 

1 



-1 



d^ / 7 --, ~ — r dz = 

(cosh q, — z) 

dz = 2Ti . (9.17) 



l + z\ 

arctan I \ I — arctan 



As result we have the following formula: 




e^-JMEE . (9.18) 

2 a 

As we see, in the case of \^ 7^ 1 this formula gives that S —>■ 00 as a —>■ 0\ Moreover 

the effective mass of this system is negative, in the sense that the acceleration (see 

( |9.15|) ) is directed to the region with a greater rest energy! 
1 2 

Notice, the sum V+V are considered here as initial approximation, in particular. 



for very small distance between the singularities (2a <^ f). However, as indicated 
above (see the end of section |^ this may be rightful, because absolute value of 
the dyon's fields D, B are infinity at the singular point. Thus we can hope that 
our iterative procedures ( |7.4| ) or ( [7.7| ) with this initial approximation can give an 



approximated solution. In any case the investigation of field configuration ( |9.1| ) may 
give leading considerations for obtaining the solution. 

Now we present the results of numerical calculations^ for the function S{a, V"^) . 
We used the dimensional system in which a = 1 and f = 1 . In this dimensional 

^Here we had used FORTRAN program with NAG hbrary (function DOIDAF for integration in two 
dimensions). 
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system the energy of for a = oo and ^ = is 2£ = 4/5/3 f» 2.4721 (see (^), 
and equality a = 1 implies that a = f in other dimensional system. According to 
( |9.13|) , ( [4.2|) we have d = db , b = (1 + a^) ^ in this dimensional system. We made 
the calculation for a = 0.00004 -h 4 and = ~ 0.95 . In order to make the best 

1 2 

convergence of the integral (both near the singular points and as r,r — oo), we 
divide the space of integration on two part. In the first part, which is bounded by 
the sphere with radius 20 -h 40 , we used the bispherical coordinates. And outside 
this sphere but inside the sphere with radius 200 -f- 1000 we used the spherical 
coordinates. The precision for obtained results of calculations is no more than 5% . 

As noted above the effective mass for the internal movement of the dyons in 
T)^^^ may be negative for very small distance between the dyons. But, of course, 
as a ^ oo this mass is positive. Thus there are an inflectional line in coordinates 
(a, V"^) on which the mass is equal to zero. We denote the point of intersection this 
line with the axes a as a . That is 

d£{a, V^) 



a 
V 



a 




(9.19) 



According our numerical calculations d ~ 4.6 ■ 10~^ and value of the energy S{d, 0) 
is near 4/3/3. It is possible that this is her exact value. This point (a, 0) in some 
sense is dual to the point (oo, 0) where the energy is £ = 2S . We shall also see this 
duality below in the context of obtained trajectories a(x°) . 

The numerical calculations give that for a < d and V = the values of the 
energy practically equal to ones given by formula ( p.l8| ). But this agreement disap- 
pears as the speed V increase. The appropriate results show in figure ^. (The values 



calculated by formula (|9.18| ) show with discontinuous line.) 
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Figure 3: Energy of D^^^ for a < a (approximation by formula ( 9. IS]) ). 
The best agreement with the numerical results is given by addition the term 



l_y2 

the condition £{d, 0) 
the energy: 



^ to expression ( |9.18| ) for the energy. Here the coefficient C is chosen from 

2 S . Thus we have the following approximate formula for 



TT 



bdVT 



2a 



2 ^ 



2^3 



1 - \/2 



(9.20) 



20 



According this formula we have the following expression: 

3 vr 6 J 



(9.21) 



In our dimensional system this formula gives a = 4.63673 ■ 10 ^ that good agree 

^_yi j , where n = 2,3,4, to 
formula ( 9.20 ) we can obtain a formula, which gives very good approximation for the 
results of numerical calculations (for a < a and < 0.95). This formula has the 
following form: 



0.0115 VI - 1.236^2 
a ^1-^2 



0.083 



0.007 



0.0002 V'^ 



Here we change the two factors in formula ( p.20|) for their numerical values in our 
dimensional system. In figure ^ is shown this good agreement between the results of 
numerical calculations and values calculated by formula ( |9.22| ). 
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Figure 4: Energy of D^^^ for a < a (approximation by formula ( 9.22D ) 



It is evident that formula ( p. 221 ) gives the same expression for a ( p. 21 



In the region of a > a the function S{a, V^) has a minimum. The appropriate 
plot for = is shown in figure ^ 
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Figure 5: Energy of D^^^ for middle values of a and V = . 

As we see, the energy falls very quick and then it increases very slow as the values 
of a increase. Such dependence is convenient to show by the plot with logarithmic 
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Figure 6: The energy in logarithmic scale of a and = 0.6 with step 0.05 . 



a = am ^ 0.045 . For very small value of deviation \a — am\ we have the following 
approximate formula: 

£ = £^ + 128 (a- a^)^ +1.651/^ . (9.23) 

Now we investigate the trajectories for various values of the energy. With the 
numerical calculations we have obtained the function V{a, £) for the region of a 
which is shown in figure |. The appropriate results are shown in figure 0. As we 

V 




Ig a 



Figure 7: The function ^(logio ^) for <S = 0.5 2.3 with step 0.2 . 
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see, here we have a periodical motion. According to formula ( p.l4| ) the period of this 
motion is given by the formula 

da , 27r , , 

and uj = — . (9.24) 



Via, £) T 

For the minimal value of the energy we calculate the circular frequency uj using 
formula ( |9.23| ). The appropriate value is 8.825. And for values of the energy from 
0.5 to 2.3 with step 0.1 we numerically calculate the frequency with help of formula 
(^1241)^- As result we obt ain the function uj{£) . The obtained dependence is 



shown in figure ||. The dual points (a, V) = (a, 0) and (oo, 0) are unstable points 
of rest with the energy 2 8 . The appropriate period is infinity and the frequency is 
zero. 




Figure 8: The function uj{£) of the periodical motion for £ < 4/3/3 . 
Now we consider the movement in the region a < d . The lines with constant 



energy for this region are shown in figure Let us consider the line in figure |T0 
that correspond to sole value of the energy. As we see in this figure, if for an initial 
instant of time the speed is zero {V = 0), then the dyons begin to approach each 
other from the points a = a/{S) and quick gather a high speed. According to the 
all formulas for the energy in the region a < a ( (|9.18| ), (|9.20| ) and (|9.22|) ) we have 
the following expression: 

a,^'-§ . (9.25) 

But at any intermediate point between a// and a, in figure the dyons may 
have two different values of the speed for one value of the energy. In principle, the 

^In fact we numerically calculated the integral from a„iin + Sai to a„iax — 5a2 ■ For the re- 
gions [amin, flmm + <5ai] , [cimax — Sa2, ttmax] we approximated the functlou £{a,V'^) by a linear 
dependence on a , and analytically calculate the appropriate integral. 
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Figure 10: Line with constant energy for a <a . 

system may jump from one branch of the constant energy hne to another. And in the 
initial approximation there are not any criterion for choosing between these branches. 
But in first approximation we take into account a radiation of the dyons. Because an 
amphtude of the radiation depend on acceleration, we can propose that this radiation 
will prevent for the jumping (with infinite acceleration). Thus we believe that the 
bidyon system in initial approximation moves from the point a — a, , V — to 
the point a,/ on the bottom (in the figure) branch of the constant energy line. At 
the point a — a,i the acceleration on the trajectory is infinity and direction of the 
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velocity is changed to opposite one. From the point a„ the dyons may move on the 
top branch ( a — oo ) or come back on the bottom one. Of course, near this point 
the radiation plays a significance role for the movement. Thus we can suppose that 
the bidyon dynamical system stop short of the point with infinite acceleration and 
come back on the branch with less speed, that is the bottom branch. In this case we 
have the periodical motion between the points a, and a„ . 

We shall numerically calculate the period for this periodical motion. But at first 
we shall have obtained some analytical result for an idealized case. If we have a 



dynamical system with energy given by formula ( |9.18| ), then we have her trajectory 
with the formula 



X 



1 - 



2a£ 
nbd 
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dd 



nbd 



arcsm 



2dS 
Tcbd 



(9.26) 



Using formula ( |9.25 ) for a, we can obtain the following trajectory: 



2 



COS —X 



irbd nh 



(9.27) 



As we see, in this case the energy is proportional to the frequency! The system 
moves from the point a = a, to the point a = and come back. In keeping with 
the conservation of full angular momentum, the dyons may not change places. At 
the point a = the speed is equal to the speed of light = 1 , for this idealized 
system. 

Now we present a dependence oj{£) that numerically calculated by formula 
( |9.22|) for the energy. As noted, this formula very good approximate the numerically 
calculated function £(a, V) . In figure |ll] we have the plot for numerically calculated 
function uj{£) (continuous line), asymptote £ = 2£ , and dependence uj{£) for 
idealized motion ( |9.27|) (discontinuous lines). 

As we see, uj ^ oo ( T ^ ) when the energy approached 2 £ from above, in 
contrast to the case when the energy approached 2 £ from below. 

We can derive also a formula that will better approximate for the numerically 
calculated function uj{£) . Let us take movement ( |9.26| ) of idealized system ( |9.18| ). 
But we shall suppose that the system stop short of the point a = and come 
back from some point O// 7^ . We shall have obtained this point from approximate 
formula ( p.2CI| ). This formula considered as equation with their associated equation 

are the system for obtaining of the point (a/,, V^) with infinity acceleration. 



d£ 



The appropriate solution is 



^V^-Kbd\f£ 



4 (£ + £ 



V 



£-2£ 
£ + £ 



(9.28) 



25 



35000 
30000 
25000 
20000 
15000 
10000 
5000 



2 . 5 



7.5 10 12.5 15 17.5 20 



Figure 11: Function uj{S) for the periodical motion in a <a. 



Using expression ( |9.26| ) we have the period for this restricted ideahzed motion in the 
following form: 
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Tihd 
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"vr 


f2a„£V 
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\ nbd J. 


.2 



(9.29) 



And for the circular frequency we have 















; arcsm 




j 




2(£ + £f 





(9.30) 



In figure |T2| is shown this function uj{£) (discontinuous curve) with the numerically 
calculated function according to formula (|9.22| ) for the energy (continuous curve). 
As we see, there are some distinction between the plots near £ = 2£ . But for 
superior value of the energy we have good agreement between the plots. For high 
values of the energy formula ( p.30| ) gives 



A£ 12V3^ 



V£ 



(9.31) 



Thus we have considered the internal movement of the bidyon for initial ap- 

1 2 

proximation, which is the sum of the dyon's fields J) + V ■ The appropriate first 
approximation is connected with the radiation of the dyons. The oscillating bidyon 
in initial approximation may give rise to some standing wave in first approximation. 
According to properties of the radiation for accelerated charge, an amplitude of this 
standing wave is zero at infinity of three-dimensional space. In this sense we can call 
it as restricted standing wave. The standing wave may change the trajectory a(x°) . 
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Figure 12: The function uj{£) for a < a with its approximation. 

Thus we shall have an energy interchange between the initial approximation and the 
standing wave. Some equilibrium oscillating field configuration may exist for this 
case. Because this problem has two-point boundary conditions, it is possible that 
the bidyon oscillating system will have some discrete set of allowable frequencies. 

Furthermore, the oscillating bidyon will interact with any plane waves, which 
are exact solutions of the model 0. The resonance interaction will intercept the 
plane waves with bidyon's frequencies. As result we shall have some non-restricted 
(in three-dimensional space) standing wave connected with the oscillating bidyon. 
Thus it is valid to say that such oscillating field configuration occupies whole three- 
dimensional space. Of course, near the singular points this configuration has a lo- 
calization region, but the associated non-restricted standing wave has some finite 
amplitude at infinity. 

We can represent the oscillating bidyon solution by Fourier series in time coor- 
dinate. Thus we have the following series for the solution in own coordinate system 
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V=Y1 2/exp(-i/c^y° 

/=-oo 



(9.32) 



where a; = — , 'D_t = 'Vf(y^,ui) , are complex column-function given by the formula 



T 



1 

2/ = — /^exp(i/^y°) di/° . 



(9.33) 



With help Lorentz transformation we can obtain from ( |9.32| ) the following moving 



27 



bidyon solution: 

oo 

V= 2^/exp(i/e) , (9.34) 

/=-oo 

where Q = k^x^^ , k^k^^ = — <J^ ■, and the fields Dj, (for the column Pj) 

are connected with the fields Dj, Bj (for the column ) by formulas of Lorentz 

transformation (|]6|a) with functions E(D, B) , H(D,B) (jgj). 



10. Effect of gravitational interaction 

In the article we had considered the distortion of a small amplitude quick- 
oscillating wave against the background of a high field, for this nonlinear electrody- 
namic model. We had shown that this distortion coincide with light beams distortion 
in gravitational field with the following metric: 

-g^^ = g^"" - O? T^^ , (10.1) 

where g^^ is the metric of basic Minkowskian space and the energy-momentum 
tensor T^^ is formulated from the background field. 

Bidyon solution (|9.34|) has the quick-oscillating wave part. As distinct from light 



wave, this wave has the rest frequency a;. But we can presume that this wave (in 
the region with small amplitude of it) will have the analogous behavior that light 
wave. 

Thus let us consider the problem for propagation of the small amplitude quick- 
oscillating wave with the rest frequency against a high field background. Here we 
shall use a method that some distinct from one, which used in the article [Q]. Let us 
take the cartesian coordinate system {x'^} with the metric g^^ = h^'^ such that 

l^ool = l^iil = 1^22! = |/i33| = 1 , /i^^ = for 7^ . (10.2) 

In this section we don't define one of two signatures of the Minkowskian metric: 

(- + ++) or (+ ). 

In this coordinate system we have the following form of equation (|2.12|) for the 
potential outside of the singular set: 

C^'^P ^ ' = . (10.3) 

Let we have a small field configuration A^{y'^) , which is periodical on with 
circular frequency u and satisfies the Maxwell equations 

C^""' ^ = , (10.4) 

dy^dy" 
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where 



This field configuration may represent the quick-oscillating part of the bidyon solu- 
tion (|9.32| ) in the region where the electromagnetic field is small and the nonlinearity 
can be neglected. With help the Lorentz transformation we can obtain also an appro- 
priate moving field configuration. But if we have its moving against the background 
of a high field, we can presume a change for parameters of this small field config- 

o 

uration. Let us take a solution of the model as the background field Af_i and a 
moving field configuration Af^{x'^) , which is related with A^{y'^) by some unknown 
transformation. We substitute their sum 



A„ 



Au +A, 



(10.5) 



into system of equation ( p.0.3| ). Taking into account that for the associated electro- 



magnetic field Ff^i, -C F^u and 
we obtain the following equation: 



dF, 



crp 



> 



d F 



p,v 



crp 



rfi A 



dxi^dx" 







(10.6) 



o 

where coefficients C^^'^'' ( |2.13|) include the background field Fpv only. 

Into some restricted region, where we can believe that the background field is 
constant, we can try to find a coordinate transformation 



(10.7) 



which transform system of equations (|10.6|) to system ( |10.4|) . We shall have this 
transformation if 



^pvap ^ QU^ QP^ (jS^a 



-g'P<r-,up 



(10. 



where g'^'' = G'^ G''^ h^^ is a metric associated with this transformation. 

With help the direct substitution we obtain that g'^^ = g'^^ , where the metric 
g'^'' is defined by formula (|10.1|) ! 

Let us write the basic relations, which connect the metric g^'' with the trans- 
formation matrixes. 



r''G':^G':, = h'^p 



=,pp 



(10.9) 



Thus we have a solution of system ( |10.6| ) into the restricted region (where 
g^^'^ fa const ) in the following form: 



A, = G^^A,iy' 



;io.io) 
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The metric g^'^ , calculated by formula ( p.0.1| ) with the background field for 



full four-dimensional space, define some curvilinear Riemann space, in general case. 
Thus, according to known theorem of Riemann geometry, we can't consider G^'^ as 
coordinate transformation tensor for full space. But we can always make the trans- 
formation to the cartesian coordinates {y^} for one point of four-dimensional space. 
As this point we take the point {x'^} , where are three-dimensional coordinates 
for some central point in a localization region of the field at the time moment 
x° . Thus for using solution ( 10.10 ) in full four-dimensional space we must take 



G'^^ = G';^{xP) in formula ([lalOD . But because x' = x'{x^) we have G^^ = C^^ix^) 
into formula ( p.0.10| ). For the coordinates {y'^} we take the following expressions: 

f = G] (x^ - x^) , G] = , x' = x'{x°) , 

y^=(k^dx^ , k,{x'') = Gl{x^) . (10.11) 



Because the solution Aj^iy'^) is periodical by ?/° , the proper time play a role of 
normalized phase for quick-oscillating wave ( |10.10| ). In this case the ordinary phase 
is uiy^ and the components of wave vector are wfc^. We have also the following 
evident relations: 



dk^ _ dku 



(10.12) 



dx" dx^ 

According to formulas (|10.9| ) and definition for k^ ( |10.11|) we have 

r''k^K = h^'' , (10.13) 
where g^" = g^'^{x'') , A;^ = k^{x'^) . 

This is dispersion relation for the quick-oscillating wave (see also 0). 

Let us adopt the proper time y^{x^) at the point x'^ as parameter of motion s 
for the localization region of the field . That is s = y^{x^) . Thus we have the 
following definition for four-velocity: 

f//^ = ^ = G^'Jx-) . (10.14) 
ds 



Definitions for ( |10.14|) , ky (|10.11 ) and formulas ( 10.9|) give also the following 
relations for the point {x^} : 

k^U^ = l , U^ = h^°g^"'ky . (10.15) 



Let us introduce the inverse tensor ^f^^ for the tensor g' 



9,^9"' = 5", ■ (10.16) 

Of course, within the limits of the introduced Riemann space the quantities are 
components of the metric tensor with inferior indexes. But for the basic Minkowskian 
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space the quantities are components of some tensor field and are components 
of appropriate inverse one. 

From ( |10.15| ) and (|10.16|) we immediately obtain for the point {x^} that 



00 



;io.i7) 



Thus we have 



— I^.dx'^dx^ 



g^^ dx^ dx" 



From ( |10.15| ) we have also 



1 



'10. 18) 



;io.i9) 



Let us obtain an equation for the trajectory x^{s) . Differentiating dispersion 
relation (^UTT^) by some coordinate x^ and using relations ( |10.12| ) and expression 
for through fc^ ( |10.15D we obtain for the point {x^} 



dxP 



dx" 







oxP ds 



. (10.20) 



Substituting ( |10.19| ) into ( |10.20| ) and using ( |10.16| ) we obtain the trajectory equa- 
tion in the following form: 



dU^' - 

WU^ 

ds 







(10.21) 



where 



■p/i 

up 



rJlS 



l5u 



/<5p 



I up 



dxP dx 



dx^ 



(10.22) 



As we see, within the limits of the introduced Riemann space (for which g^'^ = g^'^ 



9iJ.u = dpu 5 ^up = ^up ) this equation is geodesic line equation! 

Notice that here we use the primordial equation for the potential ( |10.3| ) with- 

d A 

out any gauge condition. It is evident, we can take the gauge condition g^'' " = 

for equation 



dyP 



for equation ( |10.6| ), which is transformed to the condition h^'^ 

Thus we have the Riemann space which is induced by background field. The 
relatively small quick-oscillating wave field configuration with rest frequency and 
localized near point x*(x°) moves on geodesic lines of this Riemann space. This 
looks as motion of gravitating particles. Thought the bidyon's field is high near the 
singular points, we may presume that in some cases the quick-oscillating bidyon field 
configuration will move also on these geodesic lines. We shall discuss this question 
in the following section. 
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11. Nonlinear electrodynamic world with oscillating bidyons 

In this section we shall discuss a possible world constructed from the great number of 
the bidyon-type field configurations. It is evident that in this world the bidyons and 
its possible strong coupling combinations may play the role of particles. Particles 

with full angular momentum multiples to — may be constructed from some number 
of the bidyons. It is amply clear that moving bidyon ( |9.34| ) has both particle and 
wave properties. 

Because = a <C 1 , we may build a perturbation theory with the fine struc- 
b 

ture constant a as small parameter, for some aspects of the mathematical model 
of this world. Thus we have an analogy with perturbative procedure of quantum 
electrodynamics. 

The world's field configuration is solution of the model's equations in whole four- 
dimensional space and each dyon singularity has a determined trajectory providing 
the satisfaction to the point boundary conditions. To obtain the solution for whole 
space, actually we may use the method of successive approximation only. According 
our method describing in section 0, the dyon's trajectories are corrected for each 
step of iteration with help of some integral conditions. In the present article we use 
the integral conservation law of the energy-momentum for this purpose. Thus this 
method for obtaining the trajectories is non-local in character. We have the classical 
generalized Lorentz force in initial approximation, for the case when we may consider 
the conservation law for the localization region of one dyon only. In general case the 
integral conditions may give the following effect. If we have a change of the field 
near some distant (from the trajectory) point, then we may obtain a change of this 
trajectory (without a propagation of disturbance from this distant point). This looks 
as quantum behavior when events at space separated points behave as though there 
is an instantaneous interaction between the events. According our concept, this 
pattern may be connected with the specific iterative way on which we go (in our 
calculations) to unknown world's field configuration in space-time. 

The wave vector fc^ of moving quick-oscillating bidyon ( |9.34| ) (without any ex- 
ternal given fields) satisfies the following dispersion relation: 



We have also the following relation for the vector of full energy-momentum of the 
bidyon field configuration: 



where m is the full rest energy interpreted as mass. 

However, because we have the non-restricted wave (with infinite three-dimensional 
energy integral) as a part of the bidyon configuration, how we must calculate the full 



F k, 
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energy-momentum in this case is not clear. A future consideration of the problem 
for motion of the bidyon with the associated quick-oscillating waves in given external 
field must give a value of the mass. We can denote the relation between the mass 
and the rest frequency as C = m/ui. Thus we have 

Pm = C^m • (11-3) 

To correlate between the bidyons (or multidyons) and real particles, we must 
have that C = h . To determine truth or falsity of this condition for the bidyon 
solution, we must have this solution and value of the mass. But the behavior of 
bidyon's initial approximation (|9.1| ) at very small interdyonic distances (high ener- 
gies) may suggest some optimism (see (|9.27|) , (|9.31|) ). Of course, outside of the regions 



immediately adjacent to the singular points the associated exact solution may be sig- 
nificantly different from the initial approximation. Thus this problem need further 
consideration. 

Gravitational interaction (in the sense of bimetric gravitational theory) may 
exist for the bidyon solution. The quick-oscillating part of the solution with small 
amplitude will have tendency to move on geodesic line of induced by background 
field Riemann space. We may presume that in some cases this tendency will be 
a controlling factor for motion of the bidyon, that is the bidyon will be piloted 
by its small amplitude quick-oscillating part. This case may takes place when we 
have the quick-oscillating background field such that the classical electromagnetic 
interaction for the bidyon is absent because of averaging. We must presume also that 
the averaging cancels the effects of associated resonance interaction. But induced 
metric ( |10.1| ) includes squares of the field and the appropriate effect will not wipe 
out by averaging. 

Thus we may hope that the nonlinear electrodynamic world will have a correla- 
tion with the real material world. In this case we could call an associated theory as 
Nonlinear Electrodynamic Theory of the World. Such theory may be considered as 
like Einstein's unified field theory with electromagnetic field as basic. In the context 
of this approach, all diversity of solutions and nonlinear effects for the electrody- 
namic field model must correlate with all real particles and their interactions. The 
evident practical significance of such approach consists in a possible discovery any 
new real effects with help of the theoretical investigation of the model. But whilst 
still we have more questions than answers on this arduous but feasible way. 



12. Conclusions 

Thus we have presented the initial theory for nonlinear Born-Infeld electrodynamics 
with dyonic singularities. We have proposed the method for investigation of interac- 
tion between the dyons. For a special case we have obtained the generalized Lorentz 
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force acting on the dyon. We have introduced the concept of bidyon, interpreted 
as electromagnetic model of particle with spin, and we have investigated its initial 
approximation. We have obtained the effect of gravitational interaction, which ap- 
pear within the limits of this electrodynamic model. We have also discussed possible 
correlations between the nonlinear electrodynamic world and real material world. 
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